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SUMMARY 

The effect of wide variations in chamber pressure on rocket combustor performance 

was determined for the gaseous-hydrogen - liqviid oxygen propellant system. A series of 

combustors with two different contraction ratios, 2, and 10, was used to cover a nominal 

6 7 2 

chamber pressure range of approximately 300 to 1800 psia (2. 06xl0” to 1. 24x10 N/m 
absolute) at a thrust range of 100 to 200 pounds force (444 to 888 N) . The nominal 
oxidant-fuel mixture ratio was 2.3, which is that required for maximum theoretical 
characteristic exhaust velocity. Chamber length was varied from 2 to 10 inches (5. 08 to 
25. 4 cm) at each pressure level. The chamber pressure was increased by increasing 
the propellant flow rate and by decreasing the nozzle throat area at a constant contrac- 
tion ratio. 

Performance efficiencies for the combustor with a low contraction ratio increased 
with increasing chamber pressure, while those combustors with a high contraction ratio 
decreased with increasing chamber pressure. These apparently contradictory results 
were explained by using both an experimental drop size correlation that indicates drop 
size increases with increasing chamber pressure and a vaporization model with flash 
vaporization in the nozzle. 

The analysis of the contradictory results indicated that the calculated characteristic 
exhaust velocity C* performance based on the static chamber pressure can exceed 
100 percent. This performance is possible because of the large total pressure loss that 
occurs with burning in the nozzle near the throat. 


INTRODUCTION 


As the need for higher thrust rocket vehicles increases, it becomes even more de- 
sirable to obtain a high ratio of thrust or energy release to rocket-combustion-chamber 


volume, which may be accomplished by increasing the combustor chamber pressure. 

One hindrance in the design of high-pressure combustors is the lack of controlled 
experimental performance data over a wide pressure range. It would thus be desirable 
if experimental performance characteristics were known for a single type of combustor 
for both low-pressure regions where ample experimental and analytical results are 
available (refs. 1 to 4) and high-pressure regions where only limited information is 
available (ref. 5). The purpose, therefore, of the present study is to investigate experi- 
mental performance over a wide range of chamber pressures with varying contraction 
ratios, propellant flow rates (at a constant mixture ratio), nozzle throat diameters, and 
chamber lengths. These data are compared with available analytical results. 


SYMBOLS 


A 

s/ 

C* 

Cl, C2, C3 


m 


D 


me 


mi 


D, 


cross-sectional area, in. cm^ 

chamber contraction ratio, chamber area/throat area 
characteristic exhaust velocity, + Wj) ft/sec; m/sec 

constants 

chamber diameter, in. ; cm 
diameter of liquid jet, in. ; cm 
mass median drop diameter, jam 
Dm correlated (appendix B), pm 

calculated to satisfy vaporization model (ref. 4), pm 
nozzle throat diameter, in. ; cm 
fraction of liquid oxygen vaporized 

force-mass conversion factor, 32. 17 (lb mass)(ft)/(lb force) (sec ); 
l(kg)(m)/(N)(sec^) 


Hr 


^,f 

LjO 

T,H2 

"t, H2O 


H 


T,o 


sum of sensible enthalpy and chemical energy at temperature T and stan- 
dard conditions, cal/mole; J/mole 

of fuel (H 2 ) at input temperature 

of fuel (H 2 ) at local combustion gas temperature 

of water vapor (H 2 O) at local combustion gas temperature 

of liquid oxygen (lox) at input temperature 
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4. 184x10 erg/cal; 4. 184 joule/cal 

kinetic energy of combustion gas mixture, cal/sec; J/sec 
length of straight combustion chamber, in. ; cm 


'n 

length of converging nozzle, in. ; 

cm 

‘ef 

effective length for vaporization. 

in.; 

'ef,c 

Lgj in cylindrical chamber, in. ; 

cm 

■ef,e 

Lgj at nozzle entrance, in. ; cm 


■ef,n 

Lg£ in nozzle, in. ; cm 



molecular weight, kg/mole 

input flow of fuel (H 2 ), moles/sec 

input flow of liquid oxygen, moles/sec 

oxidant-fuel mass flow ratio 

pressure lb force/in. ; N/m 

measured chamber pressure lb force/in. ; N/m 

defined by eq. (A7) 

Gas law constant 1.987 cal/(g-mole)(°K); 8.314x10^ j/(kg-mole)(°K) 

nozzle shape factor, nozzle volume/(chamber area X 

entropy at temperature T and standard conditions, j/Ckg-moleX*^) 

entropy of fuel (Hg) vapor at local combustion gas temperature 
entropy of water (HgO) vapor at local combustion gas temperature 
temperature, °K 

reduced initial liquid-oxygen temperature, dimensionless 

velocity of air, ft/sec; m/sec 

velocity of liquid jet, ft/sec; m/sec 

initial liquid-oxygen velocity, ft/sec; m/sec 

mass flow rate, lb mass/sec; kg/sec 

mass flow rate of air, lb mass/sec; kg/sec 

mass flow rate of fuel (Hg) at injector, lb mass/sec; kg/sec 

mass flow rate of liquid jet, lb mass/sec; kg/sec 
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Wo mass flow rate of oxygen at injector, lb mass/sec; kg/sec 

length from injector in chamber, in. ; cm 
lenght from nozzle entrance in nozzle, in. ; cm 
01 exponent (see eq. (A12)) 

^ defined by eq. (All) 

characteristic exhaust velocity efficiency, percent of theoretical 
y specific heat ratio 

3 3 

Pj density of liquid, lb mass/ft ; kg/m 

Subscripts: 

g total combustion gas mixture 

t nozzle throat 

1, 2, . . . refers to position as in fig. 14 

APPARATUS AND PROCEDURE 
Rocket Combustors 

The combustor dimensions, nominal chamber pressure levels, and corresponding 
propellant flow rates are given in table I. The injector, chamber, and nozzles were sep- 
arate detachable units. The combustor is illustrated in figure 1. The chambers were 
solid copper, but the nozzles were provided with internal passages for water cooling. 

The entire engine was cooled by an external water spray. 

Single-element concentric-tube injectors (see fig. 2) were used with liquid oxygen 
injected from the central tube surrounded by an annular flow of gaseous hydrogen. Dur- 
ing the program, the central tube Avas easily bent off-axis, and the tube tip had a tend- 
ency to erode. This problem was eliminated by providing the injectors with braces 
(shown in fig. 2) to position the tube and conduct heat away from the tube to the injector 
face. The inner diameter of the central tube was 0. 116 inch (0. 295 cm). The inner and 
outer diameters of the annulus were 0. 156 (0. 396 cm) and 0. 325 inch (0. 825 cm), re- 
spectively. The braces were approximately 0. 06 inch wide (0. 153 cm) and occupied 
about 24 percent of the annular flow area. The dimensions of all the injectors used for 
the program were identical. Some, however, were constructed of stainless steel and 
others of copper. 

Flat unserrated surfaces were used between the various combustor surfaces. Seal- 
ing was accomplished by means of silicone rubber gasket rings inserted in flat-faced 
grooves between the sections, as shown in figure 1. 
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Instrumentation 


Chamber pressure was measured with a strain-gage tranducer connected to the 
chamber 1 inch (0. 0254 m) downstream of the injector face. The upper frequency re- 
sponse of the chamber pressure measurement system was approximately 200 cycles per 
second (200 Hz). 

Gaseous-hydrogen flow rates were measured with a sharp-edged subsonic orifice, 
from which mass flow was calculated from measurements of upstream pressure, tem- 
perature, and pressure drop across the orifice. Liquid-oxygen mass flow rate was 
measured with a turbine -type flowmeter. The output signal of the alternating- current 
flowmeter was converted to a direct- current signal proportional to the flow rate. 

All measurements were made from strip- chart pen recorders. The accuracy in the 
characteristic exhaust velocity C* when determined with this type of system is approxi- 
mately ±1 percent. 


Operation Procedure 

Each rxm consisted of three stages: ignition, low flow, and full flow operation. Ig- 
nition was established by the hypergolic action of gaseous hydrogen with gaseous fluo- 
rine. After ignition, liquid oxygen at a reduced flow was allowed to enter the combustor. 
The gaseous -fluorine flow was then stopped, and a low pressure hydrogen- oxygen 
(Hg-Og) flame was established. The fire valves were then fully opened to permit oper- 
ation at full propellant flow. Steady-state operation during which data were taken was 
maintained for approximately 1-^ to 2^ seconds, the longer time being required at the 
lower pressures. 

Propellant feed was accomplished by means of a pressurized tank system. The 
oxygen tank and feed system were submerged in boiling liquid nitrogen exposed to atmo- 
spheric pressure. The injection temperature of the liquid oxygen was therefore at ap- 
proximately the boiling temperature of liquid nitrogen (78° K). The gaseous hydrogen 
was injected at ambient temperature (approx. 60±20° F or 290±11° K). 

RESULTS AND DISCUSSION 
General Observations 

No combustor erosion or burnout problems were encountered at chamber pressures 
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below approximately 1600 psia (1. 10x10 N/m absolute). At higher pressures, chamber 
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and nozzle erosion did occur. Damage was always greater in the combustor with a low 
contraction ratio. This difference was probably due to the fact that the chamber gas ve- 
locity varies inversely with contraction ratio. Erosion, when it did occur, took place in 
the chamber near the nozzle inlet end and in the nozzle section. The nozzle of the com- 
bustor with a low contraction ratio was severely eroded after only two runs at approxi- 
mately 2000 psi (1. SSxlo'^ N/m^ absolute) with a 2-inch-long (5. 08-cm-long) chamber. 

In contrast, the nozzle section with a high contraction ratio at the same chamber pres- 
sure survived several tests, each with the 2-inch (5. 08-cm) and 4-inch-long (10. 16-cm- 
long) chambers, and was severely eroded after only several runs with the 10-inch-long 
(25. 4-cm-long). 

The presence of the braces in the injector seemed to cause much scatter in the per- 
formance measurements of the combustor with a high contraction ratio. The braces also 
seemed to induce a low frequency instability of about 50 to 100 cycles per second (50 to 
100 Hz) in the combustor with a high contraction ratio. Therefore, the data obtained 
with the braced injector for the combustor with a high contraction ratio is not used. The 
presence of injector braces appeared to have no effect on the operation or performance 
of the combustor with a low contraction ratio. 


Performance Measurements 

The experimental results before correction for momentum pressure losses are tab- 
ulated in table II and presented graphically in figure 3. The data are grouped so that 
characteristic exhaust velocity efficiency, a percent of theoretical C*, is a function of 
measured chamber pressure with chamber length as a parameter. The efficiencies in 
table n were based on a theoretical C* calculated for gaseous hydrogen and liquid 
oxygen with the use of the method of reference 6. Most of the data points represent an 
arithmetic mean of five or more successive runs. A few points, however, represent 
single runs and are so indicated in figure 3 and in table II. 

The experimental data indicate an increase of performance efficiencies with increas- 
ing chamber pressure for the combustors with a low contraction ratio. It is observed 
that some of the efficiencies exceed 100 percent. The measurements are not thought to 
be a function of the instrumentation system, as various transducers were used to re- 
check the measurements, and the system was frequently calibrated. In contrast to the 
low- contraction- ratio data, the high- contraction- ratio data exhibited a pronounced per- 
formance loss with increasing pressure for the 2- and 4-inch-long (5. 08- and 10. 16 cm 
long) chambers. 

Momentum pressure loss . - To explain the performance above 100 percent, a mo- 
mentum pressure loss might be applied as in reference 7. If all the combustion is as- 
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sumed to occur before the nozzle entrance, a correction of +5 percent can be applied for 
a contraction ratio of 2, while virtually no correction can be applied for a contraction 
ratio of 10. This momentum pressure loss correction would still indicate performances 
above 100 percent. If, however, the extreme case of combustion near sonic velocity is 
considered, a correction of 25 percent could be applied (ref. 7). Thus for complete 
combustion, the apparent C* efficiency could be between 105 and 125 percent. It might 
be suspected, therefore, that the actual combustion is occurring somewhere between 
these two extremes. 

In view of the foregoing, the combustion model discussed more completely in appen- 
dix A was developed. The continuity, momentum, and energy equations are simulta- 
neously satisified at every location in the combustion chamber and nozzles. The results 
of the vaporization model of Priem and Heidmann (ref. 4) were used to determine the 
percent of liquid oxygen vaporized at every location in the combustion chamber and 
nozzle. Combustion in the nozzle was thus allowed, and if sufficient burning at a high 
Mach number can be obtained, large losses in total pressure will result. 

Effect of atomization model on C* efficiency . - The vaporization model (ref. 4) re- 
quires the knowledge of the mass median drop radius to calculate the percent of liquid 
oxygen vaporized. If it is assumed, as in reference 4, that the drop size is a function 
only of the liquid-oxygen jet diameter, the drop size for all the runs reported herein 
would be constant. Figure 4 shows the results of the combustion model under this as- 
sumption of a constant initial drop size of 300 microns, which gives the best fit of the 
data at low pressures. The C* efficiencies 7 ]^:^ were determined with the use of the 
calculated pressure at the experimental transducer location (1 in. from the injector 
(2. 54 cm)). It can be seen that increased with pressure for all chamber lengths 
and contraction ratios with an assumed constant drop size. Although good agreement is 
obtained for the 2-inch chamber with a contraction ratio of 2, the obtained fails to 
agree with even the trends of the other data. A constant initial drop size does not appear 
to be a good assumption for these data. 

The drop size correlation presented in appendix B predicts that the mass median 
drop diameter would increase in proportion to the square root of the chamber pressure 
at a constant oxidant-fuel ratio and injector geometry. For liquid oxygen at 78° K, gas- 
eous hydrogen at 298° K, an oxidant-fuel ratio of 2.3, and a chamber geometry as shown 
in figure 2, the relation is D = 15.45\/P~, where is in microns and is in 
psia. Figure 5 shows the results of the combustion model with the drop size correla- 
tion assumption. It is apparent that good agreement is obtained for the performance 
trends of the combustor with a large contraction ratio for chamber lenghts of 2 and 
4 inches (4. 08 and 10. 16 cm). The increasing drop size may very well be the reason for 
the decreasing performance with increasing chamber pressure. The performance trends 
for the rest of the data are, however, not accurately predicted. 
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From the preceding it would at first appear that different atomization mechanisms 
were occurring that depend on the contraction ratio and chamber length. The combustion 
gas velocity for the contraction ratio of 2 is higher than that for a ratio of 10, but the 
cold hydrogen velocity from the injector should be dominant in both cases in the region 
near the injector where atomization occurs. A change in jet atomization with increasing 
chamber length also does not seem likely. Some other phenomena appears to be domi- 
nant at high pressures and large chamber lengths. 

Flash vaporization model . - For a large chamber length, a low contraction ratio, 
and a pressure near or exceeding the critical pressure of oxygen, the liquid-oxygen 
drops remaining at the nozzle entrance may attain a temperature approaching the critical 
temperature. Reference 8 indicates that for a uniform liquid-oxygen spray the droplets 
may approach their critical temperature in the length necessary to vaporize only 5 per- 
cent of the drop mass. Upon entering the nozzle, these drops will experience a rapidly 
dropping pressure that could cause a very rapid or flash vaporization of the drops. This 
oxygen vapor then only need be properly mixed with the surrounding hydrogen gas to 
complete the combustion. These phenomena would result in increased burning near sonic 
velocity, which would increase the total pressure loss necessary for apparent effi- 
ciencies above 100 percent. 

To accommodate the aforementioned phenomena, the vaporization model (ref. 4) was 
altered by multiplying that portion of the effective length (see appendix A) which governs 

vaporization in the nozzle by e ^ where is the distance from the nozzle entrance. 
The value of a for each operating condition can be calculated from equations (All) 
and (A 12) of appendix A. These equations are an empirical correlation for of of the data 
reported herein. 

Results of flash vaporization model . - Figure 6 shows the results of the final com- 
bustion model, including the drop size correlation of appendix B and the flash vaporiza- 
tion in the nozzle. The prediction of performance trends is excellent for all conditions. 
An apparent C* efficiency, based on the static pressure at the transducer, of as high as 
107. 5 percent is obtained. The reason for these high apparent C* efficiencies can be 
seen from an examination of figure 7 in which the total pressure at the nozzle throat is 
related to the total pressure at the injector. The total pressure loss from the injector to 
the throat shovdd be reflected in tiq^ and thus the throat total pressure should be used 
in the calculation of C*. This results in curves whose values do not exceed 
100 percent, as shown in figure 8. In figure 7 it is shown that the total pressure loss 
(especially for short chambers and high pressures) greatly exceeds that predicted from 
reference 7 (negligible for j/ = 10, and 5 percent for s/ = 2). For low pressures and 
large chamber lengths, the total pressure loss approaches that of reference 7. This 
condition can be expected since the smaller drop size and large chamber length result 
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in a large percentage of the combustion occurring before the nozzle entrance and very 
little combustion occurring in the nozzle. 

The sensitivity of the calculated apparent C* efficiency to the oxidant-fuel ratio is 
shown in figure 9. The O/F effect is great, especially at high chamber pressures. 

Much of the experimental scatter may be due to the inability to hold a constant oxidant- 
fuel ratio of 2. 3. 

The need for the flash vaporization model or an alternate can be seen from fig- 
ure 10. The combustion model without flash vaporization was used to calculate the drop 
size D^. necessary to attain the measured chamber pressure. When was com- 

pared with the correlated drop size D (appendix B) and plotted against the effective 
length at the nozzle entrance, several interesting points were noted. For pressures be- 
low the oxygen critical pressure, « ^mc iwatter what the value of the effective 
length. For short effective lengths (assymptotic to » 0. 5), D^. « even for 
chamber pressures above the critical pressure. However, for large effective lengths 
and with chamber pressures above critical, * ^mc* results would indicate 

that for conditions in which the liquid-oxygen drops could not approach the critical tem- 
perature the vaporization model (ref. 4), along with the drop size correlation of appen- 
dix B, provides a good representation of the experimental results. However, when a 
sufficient amount of the liquid-oxygen drops approach the critical temperature, the model 
of reference 4 is not valid. It is quite apparent that for supercritical pressures and an 
^f at the nozzle entrance much above 0. 5, an alternate to the vaporization model of 
reference 4 must be used. 

The data points of figure 10 that are marked with an arrow represent runs in which 
the measured chamber pressure could not be attained in the calculations no matter how 
small the drop size was assumed. The calculations were thus terminated when the ef- 
fective length for vaporization at the nozzle throat exceeded 20. 


CONCLUDING REMARKS 

Performance efficiency, in terms of characteristic exhaust velocity, was measured 
for gaseous-hydrogen - liquid-oxygen rocket combustors over an extended range of 
chamber pressures, from approximately 300 to 1800 psia (2. 06X10® to 1. 24x10*^ N/m^ 
absolute) . The pressure range was covered with two series of combustors, one had a 
low contraction ratio of 2 and the other had a high contraction ratio of 10. The oxidant- 
fuel mixture ratio was held constant at approximately 2. 3 (±10 percent) which corre- 
sponds to maximum theoretical characteristic exhaust velocity. 
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The following results were obtained: 

1. The apparent efficiency of the characteristics exhaust velocity (calculated by using 
the measured chamber pressure at the injector end) of combustors with a contraction 
ratio of 2 increased with increasing chamber pressure over the entire pressure range. 

2. The apparent characteristic exhaust velocity efficiency of combustors with a con- 
traction ratio of 10 showed a decrease with increasing pressure for chamber lengths of 

2 and 4 inches (5. 08 and 10. 16 cm). Performance efficiency of a 10-inch-long (25. 4-cm- 
long) chamber increased very slightly with pressure. 

3. Calculated performance increased with increasing pressure for all combustors 
for this program, with an assumed constant initial drop size and vaporization- controlled 
combustion. These calculations were successful in predicting performance trends for 
the low but not for the high contraction ratio combustors. 

4. Calculations based on the vaporization model with drop size increasing propor- 
tional to the square root of chamber pressure predicted a performance decrease with in- 
creasing chamber pressure. These calculations successfully predicted performance 
trends for subcritical chamber pressures and for supercritical chamber pressures when 
the effective length for vaporization was low at the nozzle entrance. 

5. When the vaporization model of item 4 was modified to include flash vaporization 
and combustion in the nozzle, the performance trends for all the data were successfully 
predicted. 

6. It has been suggested previously that at high pressures a large fraction of liquid 
oxygen in a rocket combustor may attain the critical point. The results of this program 
indicate that this situation may occur when the effective length for vaporization at the 
nozzle entrance exceeds 0. 5 inch (1. 27 cm). Flash vaporization may then occur in the 
nozzle. 

7. Large losses of total pressure can occur due to combustion in the nozzle. A total 
pressure loss, between the injector and the nozzle throat, of as much as 12 percent was 
calculated with combustion in the nozzle, while a more conventional technique assuming 
isentropic expansion in the nozzle predicted only a 5-percent loss. The apparent charac- 
teristic exhaust velocity efficiency, based on the measured chamber pressure near the 
injector, can thus significantly exceed 100 percent. 

8. Combustion of liquid oxygen beyond its critical pressure can result in increased 
characteristic exhaust velocity efficiency due to flash vaporization and combustion in the 
nozzle, but this efficiency is attained at the expense of the accompanying total pressure 
loss. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 23, 1967, 

128-31-06-03-22. 
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APPENDIX A 


BURNING IN NOZZLE AND COMBUSTION ABOVE CRITICAL PRESSURE 

Combustion Model 


The following combustion model was developed to satisfy the energy, momentum, 
and continuity equations at each position in the combustor. It is assumed that the ideal 
gas equation of state may be applied and that the number of moles of combustion products 
remains constant, which follows from the assumption of the simple H 2 -O 2 reaction in 
which each mole of hydrogen consumed produces a mole of water. This assumption is 
true if combustion is complete and if the products contain no dissociation species, and 
it is valid at low oxidant-fuel mixture ratios and for completely mixed gases. 

The momentum and kinetic energy of the liquid oxygen were not considered in the 
following equations. At the injector, where the approximation is least accurate, the mo- 
mentum ratio of gas to liquid was between 4 and 20 for the data considered, and this ratio 
will increase away from the injector. Neglecting the liquid-oxygen kinetic energy was 
even less consequential than neglecting the momentum due to the high ratio of gas to liq- 
uid velocity. Only gaseous products were considered in continuity due to the small frac- 
tion of the volume occupied by the liquid. 

The energy equation as used was 


H,0 * <Nf - H, + N„U - FjH^ „ ^ KE = N,H^ , . „ 


(Al) 


The momentum equation for flow between two locations of a straight section of the com- 
bustor was 


Pl-I’2 = 


JRNf/W2T2 W^T^, 

9 I T> 

Sc^ 


The kinetic energy can be expressed as 


KE = 


JR^NjWgT^ 

2gcA^P^ 


(A2) 


(A3) 
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and for the case of sonic velocity at the throat, 


RN.y.T. 

KE. = — (A4) 
^ 2 

When the fraction of the liquid oxygen vaporized F is prescribed at each location in 
the straight section of the combustor and a pressure given at some location, the tem- 
perature and pressure can be calculated from equations (Al) and (A2) at every location. 
The values of were obtained from reference 9 and are a function of temperature 
only. 

In the converging section of the nozzle, a thermodynamic process must be defined. 
The process involves simultaneous expansion and burning. The process was approxi- 
mated by consecutive isentropic expansions and constant area combustion processes, 
since no single thermodynamic process can describe the actual process. As seen in fig- 
ure 11, the isentropic expansions from 2 to 3, 4 to 5, etc. were followed by constant 
area combustion processes from 3 to 4, 5 to 6, etc. , respectively. The equation used 
for the isentropic expansion between points 2 and 3 was 



(A5) 


No combustion occurs during this process and F is calculated at location 2. The values 
of were obtained from reference 9 and are a fimction of temperature only. 

K the combustion process is vaporization limited, the fraction of the liquid oxygen 
vaporized F can be obtained from reference 4 where the parameter controlling vapori- 
zation in the straight combustion chamber is 



where 

(P^/300)°- 

\ o, 4 Vo/100)°- 


(A6) 


(A7) 
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and is the distance from the injector (fig. 11). Within the nozzle, the governing pa- 
rameter is 


^ef, n ^ef, e 



(A8) 


Where L . is the effective length at the nozzle entrance and is L , evaluated at 

“ 61, C 

X = L , and X is the distance in the nozzle measured from the nozzle entrance. 

C/ n 

When the effective length is determined at any location, F can be determined from 
figure 12. 

When this combustion model was used, the drop sizes necessary to satisfy the ex- 
perimental pressures were consistently too high for the data for a contraction ratio of 2. 
This situation is contrary to expectations since the higher combustion gas velocity for a 
low contraction should, if anything, produce smaller drop sizes. A reexamination of 
reference 4 revealed that for heptane sprays (the only detailed information given) the 
curve of F as a function of Lef for high contraction ratios lay near the top of the error 
band, while that of low contraction ratios lay near the bottom, at least for a large 
Consequently, the upper and lower curves of figure 12 were used for contraction ratios 
of 10 and 2, respectively. These ratios are within, or are very near, the bounds for liq- 
uid oxygen given in reference 4. 


Combustion Above Critical Pressure 

When liquid-oxygen drops entering the nozzle are near the critical temperature 
(ref. 8), the suddenly decreasing gas pressure can cause a rapid or flash vaporization. 
The combustion of this oxygen vapor is probably mixing limited, but an attempt was made 
to account for the phenomena of flash vaporization by an acceleration of the vaporization 
in the nozzle. Increased combustion in the nozzle was obtained by modifying equa- 
tion (A8) to yield 



, 0. 83 Q ' 

^0. 22g0. 33 ^ 


(A9) 


Values of ot were calculated for each run so that the chamber pressure calculated from 
the combustion model was equal to the measured chamber pressure. The values of a 
were fit by the method of least squares to an equation of the form 
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/ \^3 

C«/P^\ 

(L 0.5) 
^ \300 j ® 


(AlO) 


where Cj, Cg, Cg are constants. The resulting equation was 


/3 = 0.08j^“°-^^(^^ (L f -0.5) 
\300/ ® 


(All) 


The correlation was further modified to 

a = (A12) 

V/3^ + 9 

so that exceedingly large values of ol were not obtained with large values of |3. It can 
be seen from equation (A 12) that for ^ » 3, o; = 3. 6, while for /3 « 3, of = 1. 2 /3. 

The values of Of as well as equation (A 12) are plotted in figure 14. Although there 
is some scatter in the data, equation (A12) predicts the trend of ot. 


Procedure For Use of Combustion Model 

The combustion locations referred to in this section correspond to those shown in 
figure 11. Although in figure 11 only three isentropic expansion and constant area com- 
bustion steps are shown, the calculations made in this report used seven steps. The 
three steps are shown for illustrative purposes only. More steps should represent a 
closer approximation to the actual process. The incremental lengths in the nozzle should 
be chosen smallest near the throat since this is where T and P are changing most 
rapidly. 

An approximate static pressure at position 1 (the injector) is needed first and can be 
approximated by 


P 


1 




(A 13) 


in which an apparent C* efficiency of 100 percent is assumed. The drop size for a con- 
centric tube injector can then be calculated from equation (B3) (appendix B) . If some 
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other injector is used (drop size not dependent oa pressure), the drop size may be cal- 
culated as in reference 4, and the use of equation (B3) is not necessary. 

The effective length for vaporization at each position in the chamber (eq. (A 6)) and 
nozzle (eq. (A8) or (A9)) can be calculated, and the fraction of liquid oxygen vaporized F 
at each position can be determined from figure 12. 

The energy equation (eq. (Al)) when used with equation (A4) at the nozzle throat 
(position 8) is a function only of throat temperature, which can be determined by itera- 
tion. The throat pressure is then determined by 


P 


t 



? ^t^c 


(A14) 


Equations (Al) and (A2) (eq. (A3) used for KE) will then determine T.^ and P,^, and 
equations (Al) and (A5) determine Pg and Tg. This iteration is repeated for each 
isentropic expansion and constant area burning step until the nozzle entrance conditions 
are determined (position 2). 

The temperature and pressure at any location in the straight chamber (such as posi- 
tion 1) can then be determined by the use of equations (Al) and (A2). If the combustion 
chamber is tapered, the procedure used for the nozzle must be continued for the com- 
bustion chamber. 

If the drop size is not dependent on chamber pressure (drop size determined from 
ref. 4), the calculation is completed and temperature and pressure are determined at 
each location in the engine. For a concentric tube injector, however, the new P^ must 
be used to correct the drop size and the foregoing calculation procedure repeated. The 
third repetition of the calculation loop will usually suffice, since after this Pj^ will 
change insignificantly. 

Because of the obvious length of these calculations, they were performed on the 
Lewis Research Center IBM 7090 computer. 
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APPENDIX B 


CORRELATION OF CONCENTRIC-TUBE DROP SIZE 

A study of the drop size obtained from a concentric -tube injector similar to that of 
figure 2 was made. Water was injected from the center tube and air from the annulus. 
Photographs of the resulting spray were analyzed to determine the drop size distribution 
from which the mass median drop diameter was determined. 

3 -2 2 

The variables investigated were the gas annular area (5. 65x10 to 6. 36x10 in. ; 
0. 0364 to 0. 410 cm^; 5 areas), liquid jet diameter (0. 06 and 0. 09 in. ; 0. 152 and 
0. 229 cm) and velocity (25 and 50 ft/sec or 7. 62 and 15. 2 m/sec), and relative velocity 
between the gas and the liquid jet (300 to 1000 ft/sec or 91. 2 to 304 m/sec; 4 values). 
The injector discharged into ambient air. 

The ratio of the mass median drop diameter to the liquid jet diameter Dj is 
plotted against the square root of the ratio of the gas to liquid momentum (fig. 14) . It 
can be seen that for high momenum ratios the following correlation is valid: 


fW.V^ 

1.5^/— UL 

D. Vw V 

j T a a 


Using the ideal gas equation of state and the continuity equation yields 


(Bl) 


W. 

= 1. 5 -i 
m 




4‘^a‘^a^a 


rrRPjTa 


(B2) 


To use equation (B2) for the combustion data of this report, hydrogen properties 
were substituted for air properties and liquid-oxygen properties for water properties, 
chamber pressure was used, and equation (B2) was multiplied by 0. 286. The multipli- 
cation was necessary to get the reduced drop size due to the lower surface tension and 
viscosity of liquid oxygen as compared with water. It should be noted that for a liquid- 
oxygen - hydrogen injector no significance should be attached to the absolute magnitude 
of the drop size but only to its trend with pressure. The final drop size correlation used 
was 
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= 0.485(H 
m U 





(B3) 


or with A 


^2 ~ ^ square inch (0. 4116x10 ^ m^), pj = 1. 202 grams per cubic centi 
meter, Tjj = 298° K, and O/F =2.3 


15.45 ,/pJ 


(B4) 


where P is in psia and D„ is in microns. 

^ m 

The geometric standard deviation of the water sprays when the drop size distribu- 
tions were fit to a log-normal distribution was very nearly 2. 3. 
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Chamber contraction ratio, st, 10; straight combustion chamber length, L^, 2 inches (5. 08 cm) 
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Concluded. EXPERIMENTAL RESULTS OF PERFORMANCE EFFICIENCIES TESTS 
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Braced injectors. 
^Single run. 
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Figure 8. - Calculated characteristic exhaust velocity efficiency based on 
throat total pressure against injector total pressure. 






Drop size ratio. 
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Effective length for vaporization at nozzle entrance, Lg^ g, cm 

Figure 10. - Comparison of drop size necessary for vaporization model and 
correlated drop size. 



Figure 11. - Schematic of combustion model. 


27 









NASA -Langley, 1967 33 E-3318 



